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ABSTRACT 

In this paper, we study large-scale structures from numerical simulations, paying particular attention 
to supercluster-like structures. A grid-density-contour based algorithm is adopted to locate connected 
groups. With the increase of the linking density threshold from the cosmic average density, the foam- 
like cosmic web is subsequently broken into individual supercluster-like groups and further halos. 
To be in accordance with normal EOF halos with the linking length of 0.2 in unit of the average 
separation of particles, halos in this paper are defined as groups with the linking density threshold 
pI p = 1 -\- 5 = 80, where p is the grid density, p is the average mass density of the universe. Groups 
with lower linking densities are then generally referred to as supercluster-like groups. By analyzing 
sets of cosmological simulations with varying cosmological parameters, we find that a universal mass 
function exists not only for halos but also for low-density supercluster-like groups until the linking 
density threshold decreases to 1 -I- (5 ~ 8 where the global percolation of large-scale structures occurs. 
We further show that the mass functions of different groups can be well described by the Jenkins form 
with the parameters being dependent on the linking density threshold. On the other hand, these low- 
density supercluster-like groups cannot be directly associated with the predictions from the excursion 
set theory with effective barriers obtained from dynamical collapse models, and the peak exclusion 
effect must be taken into account. Including such an effect, the mass function of groups with the 
linking density threshold 1 + 5 — 16 is in good agreements with that from the excursion set theory 
with a nearly flat effective barrier. A simplified analysis of the ellipsoidal collapse model indicates 
that the barrier for collapses along two axes to form filaments is approximately fiat in scales. Thus 
in our analyses, we define groups identified with 1 -I- 5 = 16 as filaments. We then further study the 
halo-filament conditional mass function and the filament-halo conditional mass function, and compare 
them with the predictions from the two-barrier excursion set theory. The shape statistics for filaments 
are also presented. 

Subject headings: dark matter - large-scale structure of Universe - method: statistical 



1. INTRODUCTION 

One of the key issues in cosmological studies is to 
understand the physical processes related to the struc- 
ture formation in the universe. In the cold dark mat- 
ter scenario, gravitational effects play essential roles in 
amplifying small density fiuctuations generated in the 
early universe to shape the large-scale structures seen 
today. Being directly associated with galaxies and clus- 
ters of galaxies, virialized dark matter halos have been 
widely studied theoretically and observationally. Their 
mass function, which describes statistically the forma- 
tion and evolution of dark matter halos, is shown by 
numerical simulations to follow a functional form univer- 
sally valid for a wide range of co smological models (e.g., 
iSheth et al.ll2001l:lJenkins et al.|[2001,) . Such a universal- 
ity can be largely explained in the context of halo model 
which links initial density fluctuations to nonlinear dark 
matter halos through gravitational collapse mo dels (e.g., 
IPress fc SchechtedflQTl ICoorav fc ShethllW)! ). 

Considerable efforts have been made to improve the 
spherical collapse model to include more realistic char- 
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acteristics in the modeling. It has long been realized 
that the anisotropic features contained in the initial den- 
sity fluctuations c an be magnified by nonlinear gravity 
(jZeldovich I I1970I . I1982D . It is expected that the col- 
lapse of a region first happens along the direction with 
the largest eigenvalue of the linear deformation tensor, 
thus leading to a sheet-like structure. Subsequent col- 
lapse along the direction of the second largest eigenvalue 
contracts the sheet structure to a filament. A halo can 
eventually form once further collapse occurs in the re- 
maining direction. An ellipsoidal collapse model is de- 
veloped to exten d these considerations to the nonlin- 
ear regime (e.g., ES 11971 IWhite fc Sillj IT979h . The 
peak-patch scenario further includes the external tidal 
force self-consistently into consideration and improves 
the modeling of gravitational collapse around initial 
density peaks (Bond fc Mvers,1996: Bond et all [1996) . 
ISheth et all (|200lD and lSheth fc TormenI (|2002f) incorpo- 
rate the peak-patch scenario into the excursion set ap- 
proach in an averaged way. They first obtain statisti- 
cally the averaged shape parameters of the initial tidal 
field. These averaged parameters are then used in the 
peak-patch ellipsoidal collapse model to derive the col- 
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lapse criterion. It is noticed that on average, the halo 
formation is delayed due to the anisotropy of the gravi- 
tational effects. The predicted halo mass function (MF 
hereafter) is then in good agreements with that from nu- 
merical simulations. 

Being very important in the hierarchy of large-scale 
structures, virialized dark matter halos of galaxy scale 
and above contain only ^ 40% of the total mass in the 
universe. Majority of the mass is distributed outside 
these large halos. In the language of halo model, the 
dominant fraction of the mass in the universe is contained 
in numerous small halos down to very low mass depend- 
ing on the physical properties of dark matter particles. 
These small halos present anisotropic clustering patterns 
in space, and form, together with the massive halos, cos- 
mic web structures. From the view point of the large ha- 
los, their formation and evolution are affected mainly by 
the clustering properties of the surrounding small halos 
as a whole. Thus to the zeroth order, the mass distribu- 
tion around a large halo can be described by a smooth 
component without considering the individuality of small 
halos. This approach is clearly stated in the peak-patch 
scenario (Bond & Myers 1996; Bond ct al. 1996). In the 
framework of the excursion set theory, IShen et al.l (j2006f ) 
introduce filaments and sheets to model the large-scale 
mass distribution within which virialized halos are em- 
bedded. In their analyses, filaments are treated as an 
intermediate state of the ellipsoidal collapse when the 
collapse finishes along two directions. Then these fila- 
ments represent the smoothed version of the anisotropic 
mass distribution around fully collapsed halos. 

Various approaches have been proposed to geometri- 
cally define filamentary structures in cosmological sim- 
ulations and observations. For example tessellation 
method is introduced to reconstruct the density field, 
and the edge between tess ellations naturally constitutes a 
segment of filaments fe.g..llcke fc van de Weveaerti 1987 : 



Schaap &: van de Weveaerti 120001 IPlaten et all 120071: 
Romano-Diaz fc Van de Weveaert 1 120071). The second 
order derivatives, namely, Hessian matrix, of the tidal 
field or the density field, is also widely used to classify ha- 
los, filaments, sheets and voids acc ording to the signs of 
the eigenvalues of the matrix (e.g., iHahn et al.ll2007allbl: 



Sousbie et al.' '2008'; 'Pogos van et all i^uuijt mona et ai.i 
2010a||b!; |Arag6n-Calvo ct al.l l2010D . IStoica et al.l (I2005f l 
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propose the so called Candy model for filament finding, 
in which a marked point process with a set of chosen pa- 
rameters is used to reject points at disfavored directions 
and to locate elongated filamentary segments. These 
geometrically defined filaments, however, cannot be di- 
re ctly associated wit h the excursion-set-based filaments 
in iShen et all (|2006( ) . The dynamics of long geometri- 
cally defined filaments may not be dominated by the lo- 
cal field. Therefore they may break at the saddle point 
and accrete into the two ends separately during the late 
evolution. Furthermore, many geometrical definitions of 
filaments concentrate on the features of their spatial dis- 
tribution rather than give rise to countable filamentary 
objects. 

To emphasize their dynamical structures and to com- 
pare with the results of the excursion set theory, in 
this paper, we mainly consider supercluster-like filamen- 
tary structures. We adopt a simple but natural defini- 
tion of filaments by connectivity. Specifically, we first 



obtain the density field on a set of grids from parti- 
cle positions in a simulation. Then the site percolation 
algorithm is applied to link cells together into groups 
by specifying a linking density threshold. At a high 
enough density threshold, only virialized halos are ex- 
pected to be identified. At lower thresholds, filamen- 
tary superclusters surrounding virialized halos are lo- 
cated. The global percolation of the cosmic web occurs 
when the linking density threshold reaches a lower crit- 
ical value. This is illustrated in Figure 1. From top 
left to bottom right, the linking density thresholds are 
pI p = 1 + 5 = 1,2,4,8,16,80, respectively. As we will 
discuss later, the groups identified at 1 -I- (5 = 80 corre- 
spond to virialized halos. At \ -\- 5 ~ 16, the individual 
struct ures see n in the plot are related to filaments defined 
in Sh en et all (12006). At 1 -f (S = 8, we see the global per- 
colation, and a large structure with a scale comparable 
to the size of the simulation box appears. At lower link- 
ing thresholds, the global cosmic web gets smoother. It 
should be noted however, that even for 1 -I- (5 = 1, the 
global web structure can still be seen clearly. 

The paper is organized as follows. §2 presents our 
method in detail. In §3, we analyze the mass function 
and the occupation statistics of the identified groups with 
different linking density thresholds. In §4, we compare 
our results from simulations with predictions of the ex- 
cursion set theory. Shape statistics are given in §5. §6 
contains summaries and discussions. 



2. METHOD 



To study the 
tary objects 
able numerical 



statistical properties of filamen- 
we analyze sets of publicly avail- 
imulatio ns from GIF project 
( http: / / w w w . mpa- gar ching . mpg . de . GIF) and 
project 



Virgo 

(http:/ /www. mpa-garching.mpg.de/Virgo), 
which cover a range of cosmological models and simu- 
lation parameters. In addition, we also include in our 
analyses t hree AGDM simulations kindly provided by 
Y.P.Jing (jJing fc Sutolll998D . The relevant simulation 
parameters are listed in Table 1 and Table 2. There 
are overlaps betwe en the s imulations we use and those 
analyzed in Jenkins et al.l (^2001) to derive the mass 
function of dark matter halos, and thus comparisons can 
be made directly between the two studies. The simula- 
tions we use have relatively low resolutions. However, 
they are sufficient for our purpose of study, which aims 
to investigate large filamentary objects around halos of 
galactic scale and above without concerning the details 
of individual small subhalos. 

We apply a percolation algorithm to identify con- 
nected groups. Percolation techniques have been used 
as group finders ever since the first generation of cos- 
mological simulations. The particle based percolation, 
the FOF algorithm (e.g.. lDavis et al.iil985. ). which links 
nearby particles by a given linking length, is widely used 
to locate halos. On the other hand, the site percola- 
tion (e.g..lBorner fc MolfTgslfKlvDin fc Shandarinll 19931: 
iShandarin et al.ll2004 l20l'o| r which links adjacent grids 
with density exceeding a threshold, is mainly used to ana- 
lyze large-scale structures, such as superclusters and par- 
ticularly the morphology of global percolative structures 
at the cosmic average density. In this paper, we focus 
mostly on relatively large structures between virialized 
halos and the global cosmic-average-density surface. We 




Fig. 1. — Isodensity contours with different linking density thresholds in simulation JS12. Only the largest groups (top 100 ranking in 
mass) are shown. The global cosmic web can already be seen in average density contours (1 + 5 = 1) in the top left panel. It gets sharper 
as 1 + (5 increases to 2 and 4 (top middle and top right panels, respectively). At 1 + (5 = 8 (bottom left), the cosmic web starts to break 
out, and large tree structures are seen. At 1 + (5 = 16, the web structure breaks into individual supercluster-like groups (bottom middle). 
Finally at 1 + 5 = 80, virialized halos are idenitified (bottom right). 



TABLE 1 
Simulation Parameters 



Set 


Label 


Cosmology 


^hox 


-^soften 


-^particle 


JS 


JSIO 


ACDM 


100 


0.039 


256^ 




JSll 


ACDM 


100 


0.039 


256^ 




JS12 


ACDM 


100 


0.039 


2563 


GIF 


GIF.ACDM 


ACDM 


141.3 


0.02 


256^ 




GIF.OCDM 


OCDM 


141.3 


0.03 


256^ 




GIF.tCDM 


rCDM 


84.5 


0.036 


256^ 




GIF.SCDM 


SCDM 


84.5 


0.036 


256^ 


Virgo 


VirgojVCDM 


ACDM 


239.5 


0.025 


256^ 




Virgo.OCDM 


OCDM 


239.5 


0.03 


256^ 




Virgo.rCDM 


rCDM 


239.5 


0.036 


256^ 




Virgo.SCDM 


SCDM 


239.5 


0.036 


256^ 



Note. — Lbox and Lsoften arc in unit of Mpc/i 



thus choose the latter algorithm, which has good enough 
accuracies for these large structures and can be much 
faster than the particle-based FOF operations. The site 
percolation allows us to find different types of structures 
by specifying different linking densities. These structures 
are thus enveloped by different isodensity surfaces, from 
halos at very high density regions to the global cosmic 
web at densities approaching the average density of the 
universe (see Figure 1). 
Our specific procedures are as follows. For each simu- 



lation snapshot to be analyzed, we first obtain a density 
field on a set of regular grids by CIC interpolation from 
particle positions. We set a density threshold and pick 
up only those cells with densities higher than the thresh- 
old. Then a site percolation algorithm, in which cells 
with shared surfaces are linked together, is applied to 
connect these cells into groups. For m ore details of the 
algorithm, we refer to [Newman fc Zi fi (2001 ) (see als o 
IKlvpin fc Shandarinlll993l: IShandarin et al.ii200l IMol) . 

If particles distribute uniformly within cells, the site 
percolation would be equivalent to the particle-based 
FOF with a correspondence between the linking density 
threshold for grids and the linking length for particles. 
However, particles are not spatially uniform within cells, 
thus the resolution for the site percolation depends on 
the grid size. Compact groups with size smaller than 
the grid size can be smoothed and merged artificially 
into larger groups. To test the grid effects, for different 
grid sizes, we compare halo groups identified with the 
site percolation to FOF halos with the particle linking 
length of & = 0.2. In Figure 2, we show the corresponding 
scaled mass functions vf{v) for JS12 simulation, where 
V — (5^/o'g(Af) and vf{v) — m'^[n{'m, z) / p\[d\Tim/ dhiv]. 
Here n(m, z) is the mass function for mass m and redshift 
z. The quantity (Jq{M) is the rms of the linear density 
fiuctuations at the scale corresponding to the halo mass 
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TABLE 2 

Parameters of Simulation outputs used to fit MF 



Fig. 2. — Mass functions of site-percolation groups with 1+5 = 80 
(dotted lines) and of FOF halos (circles with error bars). For the 
dotted lines, from top to bottom, the grid number is 512^, 256"^ 
and 128"^, respectively. For the case of 512^, the error bars arc 
shown. The shaded region corresponds to the range with 1 + <5 in 
between 72 and 88 (512^). 




logM (MqI 



Fig. 3. — Average overdensity of groups for 1 + <5 = 80 (top), 
1 + (5 = 16 (middle) and 1 + 5 = 8, respectively. The horizontal 
axis is the mass of groups, and the vertical axis is 1 + A =< p > /p, 
where < p > is the average density of individual groups, and p is 
the average density of the universe. 



M, and is calculated by 



d^kW'^[k,R(M)]P{k), 



(1) 



where W[k, R{M)] is the Fourier transformation of the 
top-hat window function with the characteristic scale 
R{M) = [(M/p)(3/47r)]i/3, and P{k) is the power spec- 
trum of linear density fluctuations. The value of 6c is 
taken to be 1.686. The redshift is z = 0. The symbols 
are the results for FOF groups. The dotted lines are 
for the site-percolation groups with the density thresh- 
old 1 + 5 = 80, and from top to bottom, the number of 
grids is 512'^, 256'^ and 128^, respectively. For the upper 
most dotted line, we also attach the corresponding error 
bars. The shaded region represents the range of the mass 



Label 


z 


Ho 


Ao 


F 




^box 


^particle 


TCI n 


U 


U.o 


n '7 
U. i 




l.U 


lUU 


^ n \y ^ n9 
O.U X iU 


JS12 





0.3 


0.7 


0.21 


1.0 


100 


5.0 X 10^ 


GIF^CDM 





0.3 


0.7 


0.21 


0.9 


141.3 


1.4 X lO^" 




0.5 


0.3 


0.7 


0.21 


0.9 


141.3 


1.4 X lO^" 




1 


0.3 


0.7 


0.21 


0.9 


141.3 


1.4 X lO^" 




5 


0.3 


0.7 


0.21 


0.9 


141.3 


1.4 X 10"' 


GIF.OCDM 





0.3 


0.0 


0.21 


0.85 


141.3 


1.4 X lOi" 


GIF.SCDM 





1.0 


0.0 


0.5 


0.6 


84.5 


1.0 X 10^" 


GIF.rCDM 





1.0 


0.0 


0.21 


0.6 


84.5 


1.0 X 10^° 


Virgo. ACDM 





0.3 


0.7 


0.21 


0.9 


239.5 


6.9 X 10^'' 


Note. — L^o^ 


is in 


unit of Mpch 


and 


rripart 


ale is in 


unit of M0h 



function with the density threshold from 1 -I- 5 = 72 to 
1 -I- (5 = 88 for the case of 512'^ grids. It is seen that 
at 1 -I- ^ ~ 80, the mass function of the site-percolation 
groups agrees with that of the FOF groups very well at 
the high mass end. Thus we take 1+6 = 80 as the fiducial 
threshold for virialized halos. On the other hand, the res- 
olution effect is apparent for relatively low mass halos. 
For the case of 512^ grids, halos with M < 10 ^^-^ M© 
cannot be well resolved. However, we expect the reso- 
lution effect to be weaker for larger filamentary objects 
with lower densities, which are our main concerns in the 
paper. Thus we take 512'^ as our fiducial grid number in 
the following analyses. 

For the site percolation, the density threshold is the 
critical quantity to differentiate different groups. It spec- 
ifies the overdensity level of the envelope of an identified 
group, and should be directly associated with the average 
density within the group. Figure 3 shows such a relation, 
where the vertical axis 1-t-A is, in unit of the cosmic den- 
sity of the universe, the average density within individual 
groups calculated by M/V with V the total volume of a 
group of mass M. We present the results for three sets 
of groups with the density threshold 1 -I- i5 = 80, 16 and 
8 from top to bottom, respectively. The separations of 
1 + A for the three sets of groups are clearly seen, with 
1 + A « 400, 100 and 50, respectively. It is noted that for 
halos, 1 + A « 400, rather than ^ 2 00 defined for spher- 
ical halos (e.g.. lLacev fc Coi3ll994l ). This is because the 
volume of a group calculated here is the actual volume 
occupied by the connected cells, and halos are known to 
be triaxial in shape with a typical value of ^ 0.5 for the 
long-to-short axial ratio ()Jing fc Sutol f2002') . Discussions 
in §4 show that groups with 1 + 6 = 16 can b e related to 
filam ents defined in the excursion set theory (jShen et al.l 
12006( 1. We see that they have a typical average over- 
density of ^ 100. For 1-1-5 = 8, the global percolation 
occurs, and such a cosmic web has a typical overdensity 
of ~ 50. Note that except the grid efi^ect, we do not ap- 
ply any additional smoothing for the density field in our 
analyses. 

3. MASS FUNCTION OF SITE-PERCOLATION 
GROUPS 

In this section, we analyze statistically the site- 
percolation groups, and present a generalized Jenkins 
functional form that can describe well the mass function 
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of different groups, from virialized halos to low-density 
supercluster-like groups. 

3.1. Generalized Mass Function 

It has been shown that to a very high accuracy, 
the mass function of particle-based FOF dark mat- 
ter halos follows a universal functional form, which 
is largely independe nt of cosmologica l models and 
redsh ifts (e.g., .Sheth &: TormenI 119991 : iJenkins et al.l 
l2001fl. Improved from th e original Press-Schechter form 
(jPress fc SchechtedlTOTl . two fitting formulae are widely 
used to describe such a universal mass function, namely, 
the Jenkins form and the Sheth- Tormen form. 

The Jenkins ma ss function can be written as 
penkins et al.|[200l 

i^/M = 0.5ae-l'"^+''l^ (2) 

where i/ = (5^/(Tq(M), and a, b, and c are three param- 
eters with their fitting values a = 0.315, b = 0.61, and 
c = 3.8, respectively, for halos. 

Co nsidering the ellipsoidal collapse model (jSheth et al.l 
120011 ). Sheth & Tormen (200i) derive the mass function 
from the excursion set theory, which is given by 

For dark matter halos, a « 0.615, /3 « 0.485, and a = 
0.707. 

Thus as a test for our site-percolation-based group 
analyses, we first study the mass function of dark matter 
halos identified with our algorithm. As discussed in §2, 
our percolation groups with the linking density threshold 
1+S = 80 correspond well to FOF dark matter halos with 
the linking length b = 0.2. We then define these groups 
as halos. Figure 4 shows their scaled mass function for 
all the simulation outputs listed in Table 2. The solid, 
dotted and colored lines are the simulation results for 
ACDM at z = 0, ACDM at z ^ 0, and the other cosmo- 
logical models, respectively. The heavy solid line is the 
result of the Jenkins mass function, and the dashed line 
is for the Sheth- Tormen mass function. As expected, the 
universality of the mass function for the site-percolation 
halos is clearly seen, and both the Jenkins and the Sheth- 
Tormen functional forms can fit the simulation results 
very well. 

We now turn to groups identified with lower link- 
ing density thresholds. When the density threshold de- 
creases from the halo threshold 1+S — 80, low density re- 
gions surrounding virialized halos are included in groups. 
Nearby halos can also merge into larger filamentary-like 
objects. When the density threshold approaches the av- 
erage matter density of the universe, the global cosmic 
web extending to the whole simulation box can be identi- 
fied. A natural question raised here is whether a univer- 
sal mass function also exists for these supercluster-like 
groups. To investigate this, for each simulation in Table 
2, we construct different sets of group catalogs identified 
with different linking density thresholds, from the cos- 
mic average density to the density threshold for virial- 
ized halos and to even higher thresholds. For each given 
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Fig. 4. — Mass functions of site-percolation groups with 1+5 = 80 
for simulations listed in Table 2. The solid lines are for ACDM 
model at z = 0, and the dotted linos arc for ACDM model at 
z ^ Colored lines show the results of other cosmologies. The 
heavy solid and dashed lines are Jenkins and Sheth- Tormen mass 
functions for halos, respectively. 
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Fig. 5. — The same as Fig. 4, but for 1 + (5 = 16. The heavy 
solid and dashed lines are fitted mass functions of Jenkins and 
Sheth- Tormen forms, respectively. 

threshold, we analyze and compare the mass functions of 
groups from different simulations. It is found that for a 
wide range of linking density thresholds, the universal- 
ity of the mass function remains. In Figure 5, we show 
the mass functions for groups with the density threshold 
1 + 6 = 16. The line styles for different simulations are 
the same as those in Figure 4. It is seen clearly that the 
mass function scaled with the quantity v obeys a univer- 
sal form to the level comparable to that of halos. 

To a certain extent, the universality of the mass func- 
tion for groups beyond halos can be understood qualita- 
tively in the same way as for virialized halos. For groups 
corresponding to relatively low linking density thresh- 
olds, although they have not been fully virialized, their 
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average densities are already high enough (see Figure 3) 
so that their own gravity dominates their dynamical evo- 
lution. In other words, these groups can be regarded as 
isolated structures that are in the intermediate stages 
toward forming virialized halos. Considering the process 
of gravitational collapse of an isolated region, as long as 
all the dependence on cosmology and redshift can be cast 
into the extrapolated linear density perturbations, just as 
in the spherical and ellipsoidal collapse models, the pro- 
cess can be described in a universal way. Consequently, a 
universal mass function for these groups is expected. On 
the other hand, it is also expected that the universality 
of the mass function should break down when the link- 
ing density threshold reaches a low enough level for the 
occurrence of global cosmic web structures. We find that 
the global percolation occurs at 1 + 5 ~ 8, and indeed 
the mass function for groups with that linking density 
threshold and lower does not show a universal behav- 
ior anymore. The global percolation will be discussed in 
detail in §3.3. 

The universality of mass functions for supercluster-like 
groups raises a possibility for us to find an analytical form 
for mass functions that is generalized from that of dark 
matter halos. We consider the Jenkins form of Eq. (2) 
and the Sheth-Tormen form of Eq. (3). 

As seen from Figure 4 and Figure 5, the universal be- 
havior of the mass function depends on the linking den- 
sity threshold. Thus when we fit the functional forms 
to the simulation results, we expect that the best fitted 
values for the parameters in Eq. (2) and Eq. (3) are 
functions of the linking density threshold. 

We first consider the Sheth-Tormen form [Eq. (3)]. It 
involves three param eters a, a, and /?. From the excur- 
sion set theory (e.g., iSheth fc TormenI 120021 : iShen et all 
I2OO6), the parameters a and (3 are related to the shape of 
the collapse barrier with respect to ly, and a reflects the 
overall height of the barrier. For dar k matter halos, a i s 
us ually taken to be a = 0.707 (e.g., ISheth et al.ll2001[ ). 
In lShen et al.l (|2006[ ). they extend the ellipsoidal collapse 
model to obtain the respective collapse barriers for fila- 
mentary and sheet-like objects. Their derived barriers 
for different types of objects are different only in param- 
eters a and (3. To be in accordance with their analyses, 
in our fitting here, we fix a = 0.707 and vary a and /3. In 
the next section, we will consider more general fitting to 
further discuss the relation between our results and the 
excursion set theory. 

The heavy dashed line in Figure 5 shows our best fit 
result with Eq. (3) for groups with 1 + 6 = 16. It is 
seen that the two-parameter {a and /3) Sheth-Tormen 
functional form can fit the low-mass end of the mass 
function rather well. At high mass end, however, the 
model gives a poor fit to the simulation results. This 
indicates that the simple excursion set theory cannot ap- 
ply directly to low-density groups. As we will discuss 
in the next section, this sho uld be related to th e well 
known peak-exclusion effect ()Bond fc MversI [19960 . For 
low-density groups, such effect is stronger than that for 
high-density halos, and thus the deviation between the 
Sheth-Tormen fitting and the simulations is more appar- 
ently seen in Figure 5 than that in Figure 4 for halos. 

For Jenkins form of Eq. (2), we regard it as an em- 
pirical form, and thus all the three parameters a, b and 
c are treated as free parameters in our fitting. We then 



find that Eq. (2) can fit the mass function of groups with 
different linking density thresholds very well. We further 
obtain a generalized fitting for the three parameters that 
is applicable to all the groups in consideration, from ha- 
los with the linking density 1 -|- J = 80 to low-density 
groups with 1 -I- (5 > 8. This is given by 

a = 0.5852(1 + (5)-"-i5S2 (4) 

6 ^ 0.1898(1 + (5)°-27"i (5) 
c= 1.927(1 + (5)°-i529_ (g) 

For halos with 1 + (5 = 80, we have a = 0.295, b = 0.620 
and c = 3.766 , in ag reement with the original fitting of 
iJenkins et aP (pOOl a = 0.315, b = 0.61 and c = 3.8. 
The slightly lower value of a is due to the grid effect in 
our site-percolation analyses. For l-\-S = 16, we have a = 
0.380, b = 0.401 and c = 2.944, and the corresponding 
fitting is shown by the heavy solid line in Figure 5. 

The generalized Jenkins mass function obtained here 
allows us to perform statistically the abundance analyses 
not only for halos but also for more extended low-density 
groups. Its cosmological applications will be explored in 
our future studies. 

3.2. Occupation Statistics 

Several weak lensing measurements reveal the ex- 
istence of massive dark clum ps with un u sually 
high mass-to-light ratios fe.g, lErben et alj 120001 : 
lUmetsu fc Futamasd 120001: iMahdavi et al.l 120071) . One 
proposed explanation is that those dark clumps may arise 
from the projection effect of low-density filaments with 
their elongations happening to be near the line of sight. 
Galaxies within these low density areas are thought to be 
less clustered than that in high density clusters of galax- 
ies of similar mass. However, the proper question con- 
cerned in the dark clump problem should be whether a 
filament has significantly less projected number of galax- 
ies than that of a cluster of the same projected mass. In 
other words, it is more or less the total number of galax- 
ies contained in a filament or in a cluster that matters. 

Numerical studies show that the gravitational effects 
determine domina ntly the occupation st atistics of galax- 
ies m a halo fe.g.. iKravtsov et al.l[200l . Thus the sub- 
halo occupation statistics can give important information 
on the occupation distribution of galaxies. Here we com- 
pare filament occupation distribution (FOD) of subhalos 
with halo occupation distribution (HOD) of subhalos. 
For FOD, we define groups picked up with the linking 
density threshold 1 -I- ( 5 = 16 as fil a ments , in accordance 
with the definition of ' Shen et all (|2006f ) (see §4). For 
each filament, we count the number of halos inside it. As 
discussed previously, halos with M < 10^^-^ M0 cannot 
be well resolved with our site-percolation group finder 
due to the grid effects. Thus here for occupation anal- 
yses, we use the particle FOE method with the linking 
length parameter b — 0.2 to find halos in filaments. For 
high-density virialized halos, we do not directly count the 
subhalos inside them because the simulations we used 
have limited dynamical resolutions. Instead, we adopt 
the HOD fitting result for the average nu mber of s ubhalos 
with minimum mass Mmin from iKravtsov et al.l ((2004) , 
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Fig. 6. — Occupation statistics for filaments witli 1 + 5 = 16. The 
results from simulations JSIO, JSll, and JS12 are shown. Red dots 
show the number of halos in individual filaments. Blue dots are for 
the results with subhalos added in. The red and blue solid lines 
are the average results of the red and blue dots, respectively. The 
black solid line is the subhalo HOD result in virialized halos. The 
red and blue dashed lines are the second moments of the red and 
blue dots, respectively. 

which is given by 

W=l + (^-C)''% (7) 

where M is the host halo mass, Mi/M„iin — 22, C = 
0.045 and (3, = 1.03. 

The results of the occupation distribution are pre- 
sented in Figure 6. The horizontal axis is for the mass M 
of the host filament/halo in unit of Mmin- We limit our 
analyses to AI > 10^^ Mq in accordance with the rela- 
tively low numerical resolutions of the simulations used. 
For the same reason and also concerning galactic-scale 
subhalos, we take M^in = 10^^ Mq. The FOD resuhs 
for individual filaments are shown by red dots. The red 
solid line shows the average value of N from the red dots. 
The average HOD from Eq. (7) for virialized host halos 
is shown as the black solid line. Naive comparison be- 
tween the black and red solid lines indeed leads to the 
conclusion that the number of subhalos in a high-density 
virialized halo is statistically larger than that contained 
in a low-density filament of the same mass. It should be 
noted, however, that we find halos in a low-density fila- 
ment by particle FOF group finder with b = 0.2. These 
halos can have mass well above that of the typical galac- 
tic halo, and thus are expected to further contain subha- 
los of galactic scale in them. Those subhalos cannot be 
adequately identified in our simulation analyses due to 
the limited dyna mical resolution s. On the other hand, 
in the studies of IKravtsov et al.l 1)20041 ) . they use high 
resolution simulations and their halo finder can resolve 
well subhalos and even sub-subhalos. Thus the red solid 
line and the black solid line in Figure 6 cannot be directly 
comparable. 

To make a more meaningful comparison between the 
occupation statistics of virialized halos and that of fil- 
aments, we need to add subhalos into the halos in fila- 
ments. Then for each halo found in a filament, we adopt 
Eq. (7) as the average value to randomly assign a num- 
ber of subhalos to it. The corresponding modified FOD 



results are shown in blue dots in Figure 6. The blue 
solid line is the average of the blue dots. We see that 
the blue solid line lies above the black solid line, show- 
ing that after taking into account subhalos, the average 
FOD result is actually larger by AA^ ~ -1-1 than that 
of HOD of the same mass. Such a difference may be 
understood as follows. Considering two large halos in a 
low-density filament with one of them being the largest 
halo in the filament. When the filament evolves further 
to form a virialized halo, the less massive halo is very 
likely to merge into the largest one and loses its identity, 
thus reducing the number of occupation by AN = 1. 
Therefore if there is a proportional relation between the 
subhalo FOD/HOD and galaxy FOD/HOD, the mass- 
to-light ratio for a filament is comparable and can be 
even lower than that of a cluster of the same mass, lead- 
ing to difficulties for the filament interpretation of dark 
clumps. On the other hand, although close relations be- 
tween the occupation distribution of subhalos and that 
of galaxies are expected, differences between the two can 
exist. Thus detailed analyses of galaxy occupation dis- 
tribution are further needed concerning the quantitative 
interpretation of dark clumps with filaments. 

In Figure 6, the red and blue dashed lines are the sec- 
ond moments, defined as < N{N — 1) >, of the dis- 
tributions of the red and blue dots, respectively. The 
red dashed line lies above the red solid line, showing the 
supcr-Poisson behavior for the FOD in the case without 
adding subhalos into halos in filaments. Considering sub- 
halos in halos, the blue dashed line is nearly the same as 
the blue solid line, and the distribution of the blue dots is 
consistent with the Poisson distribution. This is because 
we add in subhalos assuming the Poiss on distribution 
in acc ordance with the HOD analyses of IKravtsov et all 
(|200l . 

3.3. Critical Phenomenon 

The large-scale connected cosmic web is the most strik- 
ing feature seen in numerical simulations as well as in 
large galaxy redshift surveys. Both the power spec- 
trum of initial density perturbations and late-time non- 
linear gravitational interactions play i mportant roles in 
shaping the cosinic-we b structure fe.g.. lBond et al.iri996l : 
iShandarin et al.ll2010[) . From Fi gure 1, we can see that 
with the linking density threshold being the cosmic av- 
erage density, i.e., 1 + 6 — 1, the global web structure is 
clearly seen. As the linking density threshold increases, 
the cosmic web becomes sharper. At 1-1- (5 ~ 8, the global 
web structure starts to break out into large tree struc- 
tures with massive halos in their central regions. With 
further increased linking density thresholds, the large 
tree structures break into individual groups dominated 
mainly by their local gravity. Thus at relatively low 
linking density thresholds, the percolation groups have 
two subclasses, relatively isolated ones in low-density re- 
gions, and the large connected groups that account for 
over ~ 50% of volume occupied by all the groups (e.g., 
IShandarin et al.ll2010[ ). While the local gravity should 
still dominate the formation of isolated groups, nonlocal 
effects affect the formation of the global web structure 
significantly. Therefore there should exist a critical link- 
ing density threshold below which the universality of the 
mass function of groups breaks down due to the nonlocal 
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Fig. 7.— Merging path along 1 + (5 for halos of M ~ lO^^M© 
and M ~ 10^^ Mq, respectively. 

gravitational effects on those large tree-shaped groups. 

To understand the global percolation phenomenon, in 
Figure 7, we show the merging path with respect to the 
linking density threshold for halos of M ~ IQ'^^Mq and 
M ^ IO^^Mq, respectively. The results are for JS12 
simulation. It is seen that, for massive halos, they more 
or less stay as isolated ones until 1 + J 8. After that, 
these large groups merge into the largest structure, i.e., 
the cosmic web. This merging process is rather sharp. 
At 1 + 6^4, all the massive halos merge into the cosmic 
web, and no isolated ones are left. On the other hand, 
for relatively low mass halos, their merging path in 1 + (5 
space is extended. At 1 + (5 > 8, they gradually merge 
into larger individual groups with the decrease of the 
linking density thresholds. At 1 + S ^ 8 — 4, some of 
them merge into the cosmic web. However, there are 
still isolated halos left even at 1 + 5 = 1. The difference 
in the merging path between the massive and low mass 
halos reflects the fact that all massive halos locate at 
high density regions which eventually become parts of 
the global web. For low mass halos, while some of them 
are in high density regions, a considerable fraction of 
them are in low density void regions and can keep as 
individual groups at 1 + 5 ^ 1 . 

Figure 8 shows the scaled mass function for all the sim- 
ulations listed in Table 2. The linking density threshold 
is 1 + (5 = 8. As expected, we can see that the mass func- 
tion for different simulations behaves differently at high 
mass end, in contrary to those shown in Figure 5 with 
1 + 5 = 16. The differences are larger than the expected 
Poisson fluctuations, showing that the universality of the 
mass function does not hold anymore due to the forma- 
tion of large tree structures which later connect to form 
the global cosmic web. Thus approximately, 1+6^8 
is a critical value for the linking density threshold below 
which the global web structure starts to be apparent, 
resulting the breakdown of the universality of the mass 
function for groups. 

4. RELATION WITH THE EXCURSION SET 
THEORY 

4.1. Unconditional Mass Function 

Originated from iPress fc Schechteil (|1974[ ). the halo 
formation theory has linked virialized halos to linear 
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Fig. 8. — Mass functions for groups with 14-5 = 8. The results 
for simulations listed in Table 2 are shown. 

density fluctuations through dynamical collapse models. 
The halo mass function can then be statistically deter- 
mined by specifying a proper collapse barrier for linear 
density fluctuations smoothed over s uitable scales. To 
overcome the cloud-in-cloud problem, iBond et all (|1991[ ) 
propose the excursion set theory, which considers trajec- 
tories in i?-domain of the linear density perturbations at 
random spatial positions. Here R is the smoothing scale 
applied to smooth the linear density perturbation field. 
By relating the mass fraction contained in halos with 
mass greater than M to the volume fraction occupied by 
trajectories first crossing the specified barrier at scales 
larger than R{M), the excursion set theory can give rise 
to th e halo mass functi on with the correct " fudge factor" 
of 2 (|Bond et al.lllQQlT ). In this theory, all the nonlinear 
gravitational effects are encoded in the shape of the bar- 
rier, which in turn is determined by dynamical collapse 
models. The simple spherical collapse model leads to a 
constant barrier that is independent of the halo mass M 
or equivalently the smoothin g scale R. Taking into ac- 
count the ellipsoidal collapse. iSheth et al.l (|2001| ) obtain 
a M-dependent moving barrier for the halo formation, 
which leads to a better agreement between the derived 
mass function and that from N-body simulations. The 
moving barrier for trajectories of linear density pertur- 
bations can be written as 

= V^S,{z)[l + l3{ai^)-% (8) 

where z is the redshift, Sc is the barrier in the spher- 
ical collapse model, and a w 0.615 and /3 « 0.485. 
The corresponding mass function from the excursion set 
theory is then approxi mately determined by Eq. (3) 
(jSheth fc Tormenll2002[) . It is noted that from the el- 
lipsoidal model, we should have a = 1. However, studies 
show that a 0.7 is often required to be in agreement 
with the ma ss fun ction from numerical simulations (e.g., 
IShcth et a l.1l200l. 

[Shen et a\\ ()2006l ) further extend the ellipsoidal col- 
lapse model to consider the formation of sheet-like, 
filament-like and halo structures, defined as having col- 
lapsed along one, two and all three axes, respectively. 
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The effective barriers for the three classes of objects can 
be generally described by Eq. (8) with 

a w 0.55, /3 w -0.56 for sheet, (9) 

a « 0.28, /3 « -0.012 for filament, (10) 

and 

a « 0.61, (3 sa 0.45 for halo. (11) 

It is seen that for halos, the effective barrier increases 
monotonically with cto(M), delaying the formation of 
low-mass halos due to the nonspherical collapse. For 
sheet- like objects, it is a decreasing function of cro(Af)- 
It is for filaments that the barrier is nearly a constant in 
(To(M). 

It has been extensively shown that the excursion set 
theory with the moving barrier and the adjusted a pa- 
rameter can give rise to the mass function of virial- 
ized halos that fits the simulation results very well (e.g., 
ISheth et al.l[200ll: ISheth fc Tor^[200l (see also Fig- 
ure 4). To a certain extent, it is expected that the ex- 
cursion set theory with suitable barriers should also be 
able to model the mass function of low-density groups as 
long as their formation is dominated by their local grav- 
ity. T he sheets and filaments discussed in iShen et all 
(|2006D are among these low-density groups. The univer- 
sality of the mass function for low-density supercluster- 
like groups identified with the linking density threshold 
1 -|- (5 > 8 shown in §3 is indeed in accordance with the 
expectation of the excursion set theory. 

Here we perform quantitative studies to investigate 
if the excursion set theory can be applicable to low- 
density supercluster-like groups. Particularly, we ana- 
lyze if the filaments defined in the excursion set theory 
of S hen et all ()2006[ ) have good correspondences to the 
low-density groups found in our simulations. 

For each group catalog found in our site percolation 
analyses with the linking density threshold from 1 + 5 = 
80 for halos to 1 -I- (5 = 8 below which the universality 
of the mass function breaks down, its mass function is 
calculated and compared with the functional form of Eq. 
(3) derived from the excursion set theory. It should be 
emphasized that the parameters in Eq. (3), a, a and 
/3, are in principle, not free parameters, but determined 
by the barrier of Eq. (8) given by dynamical collapse 
models. Thus twofolds of test should be included in the 
comparison of the theory against simulations, namely the 
functional form of Eq. (3) itself, and the values of the 
parameters therein. Here we first apply Eq. (3), regard- 
ing the parameters as free parameters, to fit the mass 
function from simulations. Then the barrier with the 
fitted values of the parameters is compared with that 
expected from the ellipsoidal collapse model. In Figure 
9, we show the fitted barrier for groups with the link- 
ing density 1 + 5 — 16, labeled as 'Fitting of original 
MF". For comparison, we also show the Sheth-Tormen 
barrier for virialized halos. As expected, the barrier for 
low density groups is lower and its slope with respect 
to CTo(Af) is shallower than those of halos. However, we 
find that down to 1 + 5 = 8, no group catalogs have mass 
functions with fitted barriers that resemble the barrier of 
filaments given by Eq. (10) derive d from the ellipsoidal 
collapse model in lShen et al.l ()2006l ) . The theoretical bar- 
rier is nearly fiat in i7o(M), while the fitted barriers all 



have significant slopes. The fitted positive slope seen in 
Figure 9 reflects that for the simulation results, the sup- 
pression of the mass function at low mass end relative to 
that at high mass end is stronger than that predicted by 
the excursion set theory with the barrier from the ellip- 
soidal collapse model. The discrepancy is clearly shown 
in Figure 10, where the green solid lines are the mass 
functions with 1 + 5 — 16 for all the simulations listed in 
Table 2, and the green dashed line is the mass function 
predicted by the excursion set theory with a flat barrier. 
The disagreement between the simulation results and the 
theoretical predictions can due either to the dynamical 
collapse model that gives rise to the barrier, or to the 
excursion approach itself 

In the recent study of lRobertson et all ()2009D , they test 
the applicability of the excursion set theory for virialized 
halos against simulations. They obtain the collapse bar- 
rier directly from simulations by tracing the collapse of 
virialized regions. They conclude that while the barrier is 
consistent with that from the ellipsoidal collapse model, 
the mass function from the excursion set theory with the 
obtained barrier is not in a good agreement with that 
from simulations. This indicates the existence of some 
in trinsic shortcom ings in the excursion approach itself. 
In lMa et al.l ()2010D and lMaggiore fc Riottol (|2010aD they 
re-emphasize the importance of the non-Markovian cor- 
rections to the excursion set theory in predicting the halo 
mass function and halo bias for filters other than the 
sharp-k filter. Within the spherical collapse model, an 
analytical formulation taking into account such correc- 
tions by introducing a k, parameter is presented in Ma 
et al. (2010). Maggiore and Riotto (2010b) and Ma 
et al. (2010) also point out that the complicated halo 
formation process can be incorporated into a stochastic 
barrier to further improve the predictions of the excur- 
sion set theory. This results an additional parameter a 
to change the barrier 5c to a^^^5c and n to an in both 
the mass function and the bias for dark matter halos. 

Another problem known to the excursion set theory 
is that the mass function is derived based on the statis- 
tics of the trajectories of random points in Lagrangian 
space. On the other hand, the structure formation 
should happen mainly around peaks of the initial den- 
sity fluctuation field. Consider a peak region that even- 
tually forms a group of mass M . For a particle away 
from the central region of the peak, the gravitational 
interaction can drag the particle into the group. How- 
ever, the average linear density fluctuation at the po- 
sition of that particle obtained by applying a spheri- 
cal smoothing centered on itself over the scale corre- 
sponding to M can be lower than the collapse bar- 
rier. Thus the particle is statistically assigned to lower 
mass groups in the excursion approach. This can lead 
to over predictions of low mass groups in comparison 
with those fr om simulations (e.g.. iBond fc Mver£lll996l : 
iMonacollTggl iRobertson et al.ll2009| ). It is expected that 
such an off-center problem can have more significant ef- 
fects on low-density supercluster-like groups, such as fil- 
amentary groups, than those of high-density virialized 
halos. 

To overcome this, IBond fc MveTi (|1996[ ) propose the 
peak-patch theory for halos. In this approach, density 
peaks in the linear density fluctuation field smoothed 
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over different scales are found. Those peaks with heights 
above the collapse barrier are potential halo centers. The 
differential mass function at mass M is then related to 
the derivatives of the number of peaks with respect to the 
smoothing scale at M . Such an approach can be regarded 
as the excursion approach only on peak particles. The 
average excursion set theory on random particles in La- 
grangian space is widely used because it is believed that 
it should resemble statistically the peak excursion theory. 
However, to compare with the mass function from simu- 
lations, there is an important additional step in the peak 
excursion approach, namely, the peak exclusion, which is 
used to trim off overlapped peaks. To certain extent, it is 
expected that the peak exclusion should largely remove 
the off-center problem in the average excursion set the- 
ory discussed in the previous paragraph. To see if this is 
the case, we perform the following analyses. Instead of 
trimming out peaks in the excursion approach, we add in 
small-scale structures back to groups found in our simula- 
tions. Specifically, for each virialized halo identified with 
the linking density 1 + 5 = 80, we hierarchically increase 
1 + 5 from 80 to 240. At each level, if a new percolation 
group occurs, and it is not the most massive subgroup of 
the parent group at the previous level, it is added into the 
original group catalog as an individual group. The same 
procedures are applied for supercluster-like groups iden- 
tified with lower linking density thresholds. For example 
for groups with the linking density threshold 1 + 5 = 16, 
we hierarchically increase 1 + 5 from 16 to 80 to add in 
subgroups back into the original group catalog. We then 
analyze the mass functions of the modified group cata- 
logs and compare them with those from the excursion set 
theory of Eq. (3). In Figure 10, the mass functions for 
original halo catalogs and the modified halo catalogs for 
all the simulations listed in Table 2 are shown in black 
and blue solid lines, respectively. While the blue lines 
are systematically higher than the corresponding black 
lines at low-mass ends, the two sets of mass functions 
are not significantly different. However, for low-density 
groups, the differences are large. The green and red lines 
in Figure 10 show the mass functions for original and 
modified group catalogs with the original linking density 
threshold 1 + 5 = IQ. The red lines are considerably 
higher than the green lines at low-mass ends. Fitting the 
modified mass functions to Eq. (3), we find that for the 
original linking density threshold 1 + 5 — 16, the fitted 
barrier is consistent with a flat-shaped barrier expected 
from the ellipsoidal collapse model of Eq. (10). The a 
parameter for the amplitude of the barrier is found to be 
a ~ 0.5. This value of a < 1 may be explained by the 
stochastic barrier model proposed bv lMaggiore fc Riottol 
pOlOb) and Ma ct al. (2010). The fitted barrier is shown 
in Figure 9 labeled as " Fitting of modified MF" , and the 
corresponding mass function from Eq. (3) is shown by 
the red dashed line in Figure 10. The good agreement 
between the red dashed line and red solid lines shows 
that taking into account the peak exclusion effects (in a 
reversed way here) , the excursion set theory with t he bar- 
rier gi ven by the ellipsoidal collapse model of Shen et al.l 
(|2006| ) with adjusted a « 0.5 can describe well the mass 
function of filaments. In this sense, it is appropriate for 
us to define groups identified with the linking density 
threshold 1 -I- (5 = 16 as filaments. 
The analyses shown in this subsection demonstrate the 
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Fig. 9. — Effective barriers for the excursion set theory. The 
Sheth-Tormen barrier for halo, and the fitted barriers for the orig- 
inal and modified groups with 1 -|- 5 = 16 are shown. 




0.1 1 10 

V 



Fig. 10. — Mass functions for the original and modified groups 
of halos 1 -I- (5 = 80 and filaments (1 -|- 5 = 16) for simulations in 
Table 2 are shown. The heavy red and green dashed lines show the 
results from the excursion set theory [Eq. (3)] with fiat barriers 
fitted to the red and green solid lines, respectively. 

importance of the peak exclusion effects in the excursion 
set theory, especially when it is applied to model the mass 
function of low-density filamentary groups. Without con- 
sidering such effects, the average excursion set theory 
predicts significantly more low-mass filaments than those 
identified in simulations. 

4.2. Conditional Mass Function 

In the excursion set theory, the conditional mass func- 
tion of objects can be analyzed by invoking two sets 
of barriers that are appropriate for the two classes of 
structures in consideration. The theory has been ex- 
tensively applied to construct dark matter halo merg- 
ing trees, where /(m, zi|M, Z2), the mass fraction of the 
main halo of mass M at redshift Z2 that is contained 
in progenitor halos of mass m at redshift zi > Z2, and 
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Fig. 11. — Halo-filament conditional mass function /(m|Af). The histograms show results from simulation JS12, and the sohd linos with 
wiggles are for the results from the two-barrier excursion set theory obtained by sharp-fc Monte-Carlo simulations. 



/(M, Z2\m^ zi), the probability that a halo of m at red- 
shift zi finds itself in halos of mass M at Z2 < zi, are of- 
ten investigated. In this case, the two barriers correspond 
to the barrie rs of halo formation at zi and Z2, respec- 
tively fe.g.. iLacev fc C ole 1993; Kauffmann fc ^yh"iti 
[1991 iSomerville fcKorattiil999t iZha ng et alJ l2008D . In 
the framework proposed bv lShen et al . (2006) for the for- 
mation of different types of objects, halos, filaments and 
sheets, the conditional mass functions can also be stud- 
ied, which can potentially reveal the environmental de- 
pendence of structure formation. Here we analyze the 
halo-filament and filament-halo conditional mass func- 
tions from simulations and compare them with those pre- 
dicted from the excursion set theory. We denote m as 
the mass of halos, and M as the mass of filaments. Thus 
f{m\M) represents the halo-filament conditional mass 
function, i.e., the mass fraction of a filament of mass M 
that is contained in halos of mass m. The filament-halo 
conditional mass function is written as f{M\m), which 
gives the probability that a halo of mass m locates at 
filaments of mass M. 

As shown in §4.1, the mass function of filaments di- 
rectly from simulations with the linking density 1+S = 16 
is not cons istent with t hat fr om the nearly flat barrier 
derived by iShen et al.l (|2006l ) for filaments. A moving 
barrier with a significant slope shown as "Fitting of the 
original MF" in Figure 9 is needed. Then to test the va- 
lidity of the excursion approach itself, in our calculations 
of the conditional mass functions from the excursion set 



theory, we do not use the barriers given by IShen et al.l 
(2006) for halos and filaments. Instead, the barriers with 
the parameters a, a, and /3 obtained by fitting Eq. (3) to 
the mass function of halos with 1 -I- (5 = 80 and the mass 
function of filaments with 1 + S = 16 from simulations 
are adopt ed. For ha l os, th e fitted barrier is very close to 
that of Sh eth etaD ()2001h shown in Figure 9. For fila- 
ments, the fitted barrier is the one labeled as "Fitting of 
the original MF" in Figure 9. 

Figure 11 shows f{m\M), the halo-filament conditional 
mass function, for different values of M. The histograms 
are the results from simulations of JS12, and the solid 
lines are the sharp-k Monte-Carlo results of the excursion 
set theory with the fitted barriers for halos and filaments. 
It is seen that the excursion set theory can describe rea- 
sonably well the overall shape of f{m\M). For relatively 
small M, it overestimates f{m\M) at low m/M. This dif- 
ference should not be due to numerical resolutions of the 
simulations as we consider only halos with m > HIJ^Mq. 
Similar difference is also seen in ICole et al.l (|2008h for 
halo-halo conditional mass function, where they com- 
pare the simulation results with those from the extended 
Press-Schechter theory with the barriers from the spher- 
ical collapse model. Other studies indicate that taking 
into account the ellipsoidal collapse improves the agree- 
ment between the the excursion set theory and simula- 
tions for the halo- halo conditional mass function (e.g., 
iGiocoli et al.ll2007f) . Note that in our analyses, the solid 
lines in Figure 11 are calculated by applying the fitted 



12 



Yan & Fan 




Fig. 12. — Same as Fig. 11 but for filament-halo conditional mass function f(M\m). 



barriers from the unconditional mass functions of ha- 
los and of filaments, respectively. Thus the differences 
seen in Figure 11 should stem mainly from the excur- 
sion approach itself. The off-center problem discussed 
in §4.1 can lead to an underestimation of merging prob- 
ability and thus an overestimation of the abundance of 
low-mass objects. Although this problem is not apparent 
for /(m), the average halo mass function, it can have sig- 
nificant effects on the conditional mass function f{m\M) 
due to the relatively high-density environment. Further- 
more, the underestimation of the merging probability for 
filaments from the excursion set theory can have larger 
effects, leading to an overestimation of f{m\M) at small 
m/M especially for low-M filaments. 

In Figure 12, we present the results for the filament- 
halo conditional mass function f{M\m). The simula- 
tion results shown as histograms appear to be rather ex- 
tended. On the other hand, the results from the excur- 
sion set theory all have sharp peaks at Af/m ~ 1 . This 
should also be related to the off-center problem of the 
excursion approach. Given the halo mass m, filaments 
with M ~ m most likely contain only single halos and 
the halo centers overlap with their host filaments. Dy- 
namically, those filaments are just the extension of the 
halos therein. In such a case, the excursion approach can 
describe well the conditional mass function f{M\m) as 
expected. However, for filaments with M >> m, they 
normally have multiple halos in them, and are likely the 
merging products of progenitor filaments each with a con- 



sidered halo in it. Thus to understand the behavior of 
/(M|m), one needs to understand well the merging pro- 
cess of filaments. As we discussed earlier, the off-center 
problem in the excursion set approach is especially se- 
vere in explaining the formation of large filaments. Par- 
ticularly, in analyzing the filamcnt-halo conditional mass 
function in which the existence of halos is pre-assumed, 
we are biased to emphasize regions with relatively high 
densities. These regions have high probabilities hosting 
large filaments. Therefore even with the fitted barrier 
from the unconditional mass function of filaments, the 
excursion set theory cannot predict well f{M\m) at large 
M/m. Such a trend is also seen, albeit to a les s extent, in 
the h alo-halo conditional mass function (e.g.. iCole et al.l 
[2001 . 

Because we use the effective barriers obtained by fit- 
ting to the unconditional mass functions of halos and 
filaments from simulations, the differences between the 
excursion analyses and those from simulations for both 
f{m\M) and f{M\m) are directly related to the dif- 
ferences in the joint probability distribution f{m,M). 
In simulations, merging process to form large filaments 
generates a valley in the f{m, M) — (m, M) plot. This 
valley leads to the broad double-peak behavior for the 
filament-halo conditional mass function f{M\m) that 
cuts f{m,M) at a given m. For f{m\M) which cuts 
/(m, M) at a given M, the valley leads to a relatively 
rapid decrease of /(m|M) at small m/M. On the other 
hand, the two-barrier excursion set theory cannot give 
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rise to the valley, resulting the discrepancies seen both 
in Figure 11 and Figure 12. 

5. SHAPE STATISTICS FOR FILAMENTS 

In this section, we analyze the shape statistics of fila- 
ment groups with 1 + (5 = 16 to sec if they indeed are 
filamentary like. 

We define the shape of a group through its grid-based 
inertia tensor, which is given by 

n 

where the sum is over all grids occupied by a group 
and X is the grid central position with i,j = 1,2,3. 
Note that we do not apply density weights to grid po- 
sitions in calculating lij. Thus our measurement re- 
flects the shape of the overall spatial extension of a group, 
and the contaminations from substructures are minimal. 
By modeling approximately the spatial distribution of 
a group as a triaxial object, we can link the axial ra- 
tios of the group to the eigenvalues (Aa, Ab, Ac) of lij by 
(a, 6, c) — \/Xa, Xb, Ac, where Q >b>c are the three axes. 
Following iFranx et al.l (jl99lD . we define the triaxiality as 



— 

An object is classified as an oblate object if its T < 1/3, 
a traxial object if 1/3 < T < 2/3, and a prolate object if 
T > 2/3. 

Figure 13 presents the shape statistics. The upper two 
panels and the lower left panel are for the axial ratios c/a, 
b/a and the triaxiality T, respectively, for filament groups 
identified with the linking density threshold 1 + 6= 16 
(solid histograms) and for halos with the linking density 
threshold 1 -|- (5 = 80 (dashed histograms) . It is observed 
that filaments tend to have smaller c/a and b/a than 
those of halos. We have c/a ^ 0.45 and b/a ~ 0.6 for 
filaments, in comparison with c/a ^ 0.55 and b/a ^ 0.75 
for halos. The shape differences between filaments and 
halos are best seen in the triaxiality T statistics. For fil- 
aments, they are dominantly very prolate with T > 0.8, 
consistent with the expected configuration for filaments. 
For halos, the distribution of T is rather wide, and most 
of them are traxial in shape. The lower right panel shows 
the distribution of c/a for filaments of two mass ranges, 
M > 10" M0 (solid) and lO^^ Mq < M < lO" Mg 
(dotted), respectively. The differences are apparent with 
c/a ^ 0.35 and c/a ^ 0.5 for high and low mass fil- 
aments, respectively. From Figure 1, we can see that 
the massive filaments in the lower middle panel locate 
at the major nodes of the global cosmic web (see for ex- 
ample, the upper right panel of Figure 1). They usually 
contain multiple large halos, and their spatial extension 
reflects the spatial orientation of the global web struc- 
ture. Dynamically, along the direction of the major axis 
of these massive filaments, collapse has not happened 
yet. Thus to a large extent, these filaments should corre- 
spond closely to the two-axis collapsed filaments defined 
from linear density fiuctuations. For relatively low mass 
filaments, most of them are extended structures of indi- 
vidual halos of similar mass. Therefore their shape dis- 
tribution should be in accordance with that of halos. On 



the other hand, because they include extended and dy- 
namically unrelaxed regions, they tend to be somewhat 
more filamentary-like than halos therein. 

6. SUMMARY AND DISCUSSION 

Applying a grid-based site percolation method to nu- 
merical simulations, we study groups identified with dif- 
ferent linking density thresholds 1 + S. Groups with 
1 + S = 80 correspond well to FOF dark matter ha- 
los. Lowering 1 + 6 allows us to find supercluster-like 
groups beyond virialized dark matter halos. As the link- 
ing density threshold approaches the average density of 
the universe, the global cosmic web structure can be nat- 
urally found. In the studies presented in this paper, we 
focus on supercluster-like groups, which are expected to 
be dynamically bound, although not virialized yet. These 
groups provide immediate environments to dark matter 
halos therein. Therefore understanding their properties 
is an important step towards understanding the environ- 
mental effects on the formation and evolution of galaxies. 

Our analyses reveal that similar to dark matter halos, 
the mass functions of supercluster-like groups for differ- 
ent simulations listed in Table 2 also follow a universal 
behavior. This universality is consistent with the consid- 
eration that these groups are gravitationally bound sys- 
tems, and form mainly through their own gravitational 
interactions. In other words, the universality found for 
supercluster-like groups and that for dark matter halos 
should arise from the same origin. We further find that 
the Jenkins functional form can describe well the mass 
functions for not only halos, but also supercluster-like 
groups. An extended Jenkins mass function applicable 
to both halos and supercluster groups is then explicitly 
presented, in which the parameters a, b, and c depend 
on the linking density threshold 1 + 6. As expected, the 
universality of the mass functions breaks down for groups 
with the linking density 1 + 6 < 8 where the global web 
structures occur. 

We also compare the mass functions from simulations 
with those from the excursion set theory with effective 
barriers derived from the ellipsoidal collapse model. For 
halos with 1 + 6 = 80, consistent with other studies, the 
two agree very well with the parameter a adjusted to 
be a = 0.707 for the moving barrier. For supercluster- 
like objects, the ellipsoidal collapse model gives rise to a 
nearly flat barrier for filaments defined as two-axis col- 
lapse objects (e.g. Shen ct al. 200ii). However, incorpo- 
rating this barrier into the excursion set theory predicts 
a mass function that cannot fit to any mass function 
of supercluster-like groups identified in simulations with 
the linking density threshold 1 -I- i5 < 80. The off-center 
problem in the excursion set theory leads to a signifi- 
cant over prediction for the mass function at low mass 
end. Taking into account this problem in the comparison, 
we find that the mass function of the groups identified 
with 1 -I- (5 = 16 is in good agreement with that from 
the excursion set theory for two-axis collapse filaments. 
Defining these groups as filaments, we further study the 
halo-filament and filament-halo conditional mass func- 
tions. Deviations from the predictions of the two-barrier 
excursion set theory are seen, which are especially signif- 
icant for filament-halo conditional mass function. 

The studies carried out in this paper can have impor- 
tant cosmological applications. The universality of the 
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Fig. 13. — Shape statistics for halos and filaments with mass largi 
and lower right panels, the dotted histograms are for halos, and solii 

a possibility for us to probe cosmologies with superclus- 
ter abundances. It can also be applied to model sta- 
tistically how the projection effects affect clusters' weak- 
lensing signals. In the very recent paper bv lMurphv et al.l 
(|2010f l. they identify filamentary galaxy groups from the 
2dfGRS survey using galaxy FOF method, and compare 
the properties of the groups with those of mock sur- 
veys constructed from numerical simulations. This study 
shows that it is becoming feasible observationally to ana- 
lyze filamentary galaxy groups statistically, and they are 
in turn can be used as cosmological probes. Physically, 
we expect that these filamentary galaxy groups should 
be closely associated with the supercluster-like dark mat- 
ter groups in our studies. To quantitatively understand 
the relation between the two, detailed FOD modeling for 
galaxies in supercluster-like dark matter groups is nec- 
essary. We discuss the FOD for subhalos in §3.2. For 
galaxies, however, the FOD can be much more compli- 
cated, and thorough investigations are highly desired. 

It is further noted that analyses on real galaxy groups 
can only be done in redshift space. Redshift distortions 
from peculiar velocities of galaxies can affect group iden- 
tifications, and further their mass functions and shape 
statistics considerably. For supercluster-like groups, 
their ambient member galaxies tend to be in the stage of 



3r than IO^^Mq for simulation JS12. For the upper left, upper right 
i histograms are for filaments. In the lower right panel, the solid and 

< M)helft'fiMofali^paiid</itlbrus their distribution suffers oblate 
distortions in redshift space. On the other hand, for 
their virialized inner regions, the distortion can generate 
finger-of-God structures in redshift space. The detailed 
impacts of redshift distortions on entire supercluster-like 
groups will be explored in our future studies. 
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